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Abstract 

This paper is devoted to the eigenvalue problem for the quantum Gaudin system. 
We prove the universal correspondence between eigenvalues of Gaudin Hamiltoni- 
ans and the so-called G-opers without monodromy in general gl^ case modulo a 
hypothesys on the analytic properties of the solution of a KZ-type equation. Firstly 
we explore the quantum analog of the characteristic polynomial which is a differ- 
ential operator in a variable u with the coefficients in {7(gl„)®^. We will call it 
"universal G-oper". It is constructed by the formula " Det" {L(u) — d^) where L{u) 
is the quantum Lax operator for the Gaudin model and "Det" is appropriate defi- 
nition of the determinant. The coefficients of this differential operator are quantum 
Gaudin Hamiltonians obtained by one of the authors (D.T. hep-th/0404153 ). We 
establish the correspondence between eigenvalues and G-opers as follows: taking 
eigen-values of the Gaudin's hamiltonians on the joint eigen-vector in the tensor 
product of finite-dimensional representation of gl^ and substituting them into the 
universal G-oper we obtain the scalar differential operator (scalar G-oper) which 
conjecturally does not have monodromy. We strongly believe that our quantization 
of the Gaudin model coincides with quantization obtained from the center of univer- 
sal enveloping algebra on the critical level and that our scalar G-oper coincides with 
the G-oper obtained by the geometric Langlands correspondence, hence it provides 
very simple and explicit map (Langlands correspondence) from Hitchin D-modules 
to G-opers in the case of rational base curves. It seems to be easy to generalize the 
constructions to the case of other semisimple Lie algebras and models like XYZ. 
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1 Introduction 
1.1 Gaudin system 

Gaudin model was introduced in |Gaudinj (see section 13.2.2) as a limit of the famous 
XXX-Heisenberg model, which describes interaction of spins in a one dimensional chain. 
Gaudin introduced his model for arbitrary semisimple Lie algebra, while original physi- 
cal motivation was restricted only to the s[{2) case. Recently spin chains with algebras 
other then s[(2) has found unexpected applications in gauge theories (see |Belitsky04| 
for recent comprehensive survey). Gaudin model also attracts much attention in mathe- 
matics due to its relation to geometric Langlands correspondence |Frenkel95l IFrenkel04] , 
Knizhnik-Zamolodchikov and isomonodromy deformation theory |L097j . Hitchin system 
!Enriquez95! , |Nekrasov95j ; geometry of polygons |KM9^ IFM03j . 

Let us recall the Gaudin construction and some notations which will be used in the 
main text. Let us denote by $ the following n x n-matrix with coefficients in we put 
element Cij on ij-th place of the matrix. We also consider the direct sum 0t„ © ••• © st„ 
and denote by $j the matrix defined as above but with the elements from the i-th copy 
of g[„ in [„©...© Let us introduce the Lax operator for the Gaudin model: 

L{u) =Y.^ (1) 

where a,- G C. 
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^ - are called quadratic Gaudin hamiltonians 

[Zk — Zj) 



To define the classical Gaudin model one needs to describe the classical phase space of 
the system. This is taken to be Oi x ... x O^, where Oi are coadjoint orbits with Kirillov's 
symplectic form. Functions TrL^{z) can be restricted from gl* to coadjoint orbits and 
they give rise to completely integrable hamiltonian system - the Gaudin system. 

In the quantum case one considers the same formulas for the Lax operator and 
quadratic hamiltonians only changing the generators of the Poisson algebra S'(g[„)®^ 
by the same generators of the universal enveloping algebra [/(gt^)*^^. It is well-known 
that Hk commute on classical and quantum levels. On classical level one can see from 
r-matrix technique that TrL^{z) and TrL^{w) Poisson commute for arbitrary k,p,z,w. 
So taking TrL^{u) for k = 1, ...,n and taking coefficients of these expressions at ^^_],,^m 
one obtains the set of algebraically independent hamiltonians. This is not true in general 
on the quantum level, for example one obtains: [Tr{L'^{z)),Tr{L'^{u))] ^ CRT04]. 

Mathematically speaking Gaudin hamiltonians generate a large commutative subalge- 
bra in t/(gl„)®^ which can be easily completed to maximal commutative subalgebra. 

To define the quantum Gaudin model it is also necessary to describe the Hilbert space 
of the model. This is taken to be (in accordance to Kirillov's orbit method ideology) the 
tensor product of some representations of g[®^: Vi ... ® Vat. Hamiltonians act on 
this Hilbert space and the quantum problem consists in finding their spectrum, matrix 
elements and so on. 

The integrability on the quantum level is a corollary of the result of Feigin and Frenkel 
|FF92j . They proved the existence of a large center for universal enveloping algebra 
of the Kac-Moody algebra on the critical level. However their result does not provide 
explicit formulas for generators of the center on critical level (due to problems with normal 
ordering), nor the explicit expressions for quantum hamiltonians of the Gaudin model. 

The explicit formula for quantum commuting hamiltonians for the Gaudin model was 
found in jTalalaevOlj : 



Qhiu) = Tri,...,„A„(Li(M) - 9„)(L2(m) - 9„) . . . iUiu) - 0^)1 (3) 



where L{u) is the Lax operator for the Gaudin model, i.e. it lies in Matnxn ® f^(0l„)®^, 
An is the antisymmetrizer in (C")®", the notation Li implies that matrix components of 
L{u) lie in i-th component of the tensor product of Matnxn- 

This result was obtained using the Bethe commutative subalgebra in Yangian algebra 
F(g[^). This subalgebra was invented in Faddeev's school, it is mentioned in [SklyaninOSl 
formula 3.16 and under the name "Bethe subalgebra" appeared in |JNU95j (where it was 
attributed to Kulish,Sklyanin,Kirillov,Reshetikhin). 

Remark 1 Let us mention that the valued matrices $ has been also considered in 
iKirillovOOl and interesting conjectures about them were formulated. 

One can take the other phase space for the same hamiltonians: the direct product of 
T*GLn, then the obtained model will be the Hitchin model on the curve CP^ with N/2 
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double-points, mixing these two phase spaces together and using reduction one can obtain 
the trigonometric Calogero model (see |Nekrasov95t ICT03] ). 



1.2 Main results 

Let us introduce differential operators with values in ?7(0[„)®^ 

QDi{u) = Tn_^An{Li{u) - du){L2{u) - 9„) . . . {U{u) - d^) (4) 
which in the case i = n 

QD4u) = ''Def'{L{u)-du) 

we will call the universal G-oper hy the reasons explained below. Here the usual notations 
of formula Q are used. One can easily see that TrAnXi...Xn = Det{X) in the case of 
a matrix X with commutative entries. The expression above is a kind of " quantization" 
for the characteristic polynomial: Det{L{u) — A). In the theory of integrable systems this 
can also be called "quantum spectral curve". The fact that A has been substituted by the 
du is related to the fact that in separated variables the Poisson bracket is {u, A} = 1. 

Remark 2 One can see that for 2x2 matrix 

the deformed determinant is just the symmetrized ordinary determinant 

TrA2{Xi){X2) = ^((TrX) - TrX^) = ^{ad + da - be - eb) 

(pay attention that in our case the symmetrization is more complicated due to the de- 
pendence on u and du)- 

The first important theorem (see section |2]) consists in the following: 

Result 1 All eoeffieients of the differential operator QDn{u) are quantum hamiltonians 
of the Gaudin system, and the following formula holds: 

TiiTi 1 ) 

QDniu) = QIn{u)-nQh^^{u)du+ ^ ^ ' QIn^2{u)dl-... 

n 

= Y.^-irGl,QI,Uu)dl (5) 
The second result is: 

Result 2 Substitution of coefficients QDn{u) by their images in the tensor product of 
finite- dimensional representations of GL{n) gives matrix-valued differential operator with 
singularities at points u = Zi, which however does not have monodromy modulo a hy- 
pothesys formulated below on the analytic properties of the KZ-type equation. 
As a corollary one obtains the following: 
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Result 3 Substitution of coefficients QDn{u) by their eigen-values on their joint eigen- 
vector gives a scalar differential operator (scalar G-oper) without monodromy. The 
construction and the conjecture above were anticipated in the hterature. The origin of 
such claims goes back to R. Baxter, M. Gaudin, further E. Sklyanin has related this to 
the separation of variables. Recently some close results were obtained by E. Frenkel, 
E. Mukhin, V. Tarasov and A. Varchenko |FMTVj . Our main source of knowledge and 
insparation on the subject is Frenkel95] and the conjecture in this form is motivated by 



this paper. 

The importance of this result is based on the following fact: in this way one can easily 
obtain the Bethe ansatz equation i.e. to find the spectrum of the Gaudin model. We will 
do it in the next paper. The idea is explained in |Frenkel95] where the scalar G-oper was 
described explicitly only for q12, gl^. 

Motivated by the analogy with Frenkel's paper we strongly believe that result El is 
true in both directions: the scalar G-oper does not have monodromy if and only if its 
coefficients are eigen-values for the Gaudin hamiltonians. 

These results are strongly connected with the geometric Langlands correspondence of 
Beilinson and Drinfeld, and Sklyanin's separation of variables theory, as it is explained in 
|Frenkel95] . We recall these connections in section |3] 

Remark 3 In the Yangian case (it is related to the XXX spin-chain) we also introduce 
the concept of the universal G-oper: 



jjQYang ^ T T A^{e-^^'-T^{u) - 1 ) (c" ^^"Ta (m) - 1) . . . (e-^^«T„(u) - 1) (6) 

This expression corresponds also to the "quantization" of the characteristic polyno- 
mial, but here A is quantized to e This is in full agreement with the Poisson bracket 

between separated variables u and A. 

One has an obvious known fact: 

Result 4 All coefficients of the difference operator UG^°''^^{u) are quantum commuting 
hamiltonians of the XXX spin chain system. 

By analogy with the Gaudin case we state the following: 

Conjecture 1 Scalar difference operator obtained from UG^"'^^{u) by substituting its 
coefficients by their eigen-values on their joint eigen-vector in finite- dimensional repre- 
sentation of Yangian has only univalued solutions. 

The scalar discrete G-opers were introduced recently by Mukhin and Varchenko see 
refrence 2 in |FMTVj . we believe that our scalar Yangian G-oper coincides with their. 



1.3 The plan 

The proof of result [T] is based on the ideas of |Talalaev04] and use of Yangian technique, 
especially Bethe subalgebra. 

The proof of monodromy properties of the G-oper takes the following lines: 
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• We prove the formula 



A^{L^{u) - a„) . . . (L,(m) - d^)A^ = n!A„(Li(M) - 5„) . . . (L,(m) - d^) (7) 
With the help of a solution of 

(9„ - L{u))~S{u) = (8) 

we define 

^>{u) = An{vi® ...®vn-i®S{u)) (9) 
which provides a solution for 

QDn^{u) = (10) 

for Vi - generic vectors in C". Moreover the corresponding linear differential equation 
has regular singularities and hence the monodromy is controlled by residues of L{u) 
considered as elements of Matnxn i^) f/(gl„)®^. 

• Unfortunetely our case is exactly the resonance case and we only state a hypothesys 
that the linear equation has no monodromy which we have verified in some particular 
cases. 

• The final observation consists in restricting the global solution of the generalized Co- 
oper to the eigen vector of Gaudin hamiltonians obtaining the scalar G-oper without 
monodromy. 

Acknowledgements. We are indebted to A. Gorsky, S. Kharchev, M. Olshanetsky 
and L. Rybnikov for the discussions and to A. Varchenko for pointing us the gap in the 
preliminary version of the paper. The work of both authors has been partially supported 
by the RFBR grant 04-01-00702. 



2 Universal G-oper and quantum commuting hamil- 
tonians 

Let us recall standard Yangian notations. Y{gi^) is generated by elements t^j'' subject to 
specific relations on the generating function T{u, h) with values in Y{qI^) EndiU^) 

T{u, h) = j^Eij ® ti,{u, h), u,{u, h) = 6,, + ^^y«-'=h^ 

i,j k 

where Eij are matrices with 1 on the i, j-th place. The relations on this generating 
function involve the Yang it!-matrix 

Pi 

R{u, h)^l - -J^^ij^ 
5 



and take the RTT form: 

R{z - M, h)Ti(;z, h)T2(u, h) = \\)Ti{z, \\)R{z - u, h) 
in Enrf(C")®2 ® Y{q[^)[z,u, h], where 

Ti(z, h) = ^ E,j ® 1 ® tij{z, h), T2(m, h) = ^ 1 ® tij{u, h). 



There is a known reahzation of Y{qI^) in f/(0t„) i.e. a homomorphism pi : 

T{u, h) = l + -y"Eij®eij = l + -, (11) 



where are generators of The following expression also provides a realization of the 
Yangian in U^qX^®^ for a given /c-tuple of complex numbers a = (zi, . . . , z^^') 

T"{u, h) = T1(m - h)T\u -Z2,h)... T''{u - zn, h), 

where T'(m — h) is the realization given by ()lip with lying in the l-th component of 

?7(g[„)®^. Let pa be the corresponding homomorphism pa : ^(fl^n) — > f^lfl^n)*^^- 

The Bethe suhalgehra in y(g[„) is generated by coefficients of the following expressions 

r,(u, h) =TrA„Ti(u, h)T2(u-h, h)...T,(M-h(z-l),h) i = l,...n. 

One has 

[ri(u, h),r,(w, h)] =0 
In |lalalaevn4j it was considered formal expressions 

i 

S^ = J2^-lyCiT^~, (12) 
j=0 

where tq = TrAnl = n\ and it was proved that 

Si{u,h) = h'QIi{u) + 0{h'+') 
where Qli are quantum Gaudin hamiltonians given by the formula 



QIi{u) = Tri_nAn{Li{u) - du){L2{u) - a„) . . . {Li{u) - 9„)1 



Theorem 1 All coefficients of the differential operator QDn{u) given by (0) are quantum 
hamiltonians of the Gaudin system, and the following formula holds: 

QDr,{u) = QIniu)-nQIn-l{u)du + ^^^QIn-2iu)dl-... 
n 

= Y.{-irCl,Qh^,{u)dl (13) 
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Proof Let us recall that the expression QDi can be obtained as a coefficient at h* (which 
is in fact the lowest non trivial coefficient) of the following expression 

TrA„(e-h^"Ti(n) - l)(e-h'9«T2(n) - 1) . . . {e-^^^T,{u) - 1) 

i 

= ^r,(?z-h,h)(-l)^-^C/e--^'h^" (14) 

j=0 

Then expanding the exponential series and collecting terms with we obtain 

J2T,iu-h,h)i-ir~^ci^-^^dt (15) 

i=o 



Next we observe the following: 

■k ^ 



i! '''"'(j-m)! 



where the coefficients -Bfc,m solve the system of recurrent relations 

— l,m— 1 

and do not depend on j. Hence 

\k 



" k\ ' ' ^ '''''' "in- m)\k\ 
Then the formula (fT^ can be rewritten as follows: 



E f - h> h)(-l)"-^'=C- ) 

m=l ^ '' ' \j=0 J 

Further one needs to note that 

n 

Y^r,[u-\\^){-\r-^Cl-J^ (16) 

= TrA„e-h^«ri(«) . . .e-h^«T„(«)(e-h^«T„+i(«) - 1) . . . (e-h^«r„(«) - 1)1 

and has the following expansion on h 

h"-™Q4_^ + 0(h"-™+i) 

Then collecting all terms of (|15|) of lowest degree in h namely at h" and using the fact 
that = 1 we obtain the theorem 



"^We apologize for this too technical proof of this result, we are sure that there exists much more simple 
and conceptual demonstration. 
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3 Solution and monodromy of G-oper 
3.1 Antisymmetrization 

We construct a global solution of ()13p in this section recalling that this differential operator 
is obtained be taking the trace of the lowest coefficients in h of a more simple expression 
which provides a factorization: 

X{u) = A4e-^9-T,{u)-l){e-^^-T2{u)-l)...{e-^^-T4u)-l) 

= h"A„(Li(^i) - du)iL2iu) -du)... (LM - du) + 0(h"+i) (17) 

Proposition 1 

A^ie-^^-T^{u) - l){e-^^-Uu) - 1) . . . {e'^^^TM - l)An = 

n!A„(e-h^"Ti(«) - l)(e-h^«T2(«) - 1) . . . {e'^^-TM - 1) (18) 



The proof of the proposition follows immediately from the lemmas below. As a corollary 
we obtain the similar proposition for the case of the Gaudin Hamiltonians considering the 
coefficient at of expression above. 



Corollary 1 



An{Li{u) - du){L2{u) - du) ... {L„{u) - (9„)A„ = 

n\ArXL^{u) - du){L2{u) - du) . . . {Ln{u) - du) (19) 



Lemma 1 // X(... ® a ® h ® ...) = — X(... 6 a ® ...), then n\X = XAn, where 
(... ® a (g) 6 ® ...) is an arbitrary vector m (C")®", X G ^® Mat®xn, ^ 

is an arbitrary as- 
sociative algebra, Mat is the algebra of matrices and An € Mat®^ is the antisymmetrizer 

Proof Arbitrary permutation can be obtained as the product of the transpositions ex- 
changing only the neighbors i and i + 1, so the conditions of the lemma above imply 

X{vi 0...0Vn) = (-l)'^"(")X(t;,(i) (g> ... O V^^n)). 

Hence 

XAn{Vi O ... ® Vn) =X{J2 (-l)'^"^"H^<x(l) O ••• ® ^a(n)) = n\X{vi (g) ... ® V^) 

adSn 

□ 

Lemma 2 

A„(e-h'9„j^i(M) - l)(e-h^«T2(M) - 1) . . . {e-^^^Tn{u) - 1)(... ®a®h® ...) 
= -A„(e-h^«Ti(M) - l){e-^^^T2{u) -I)... {e-^^^Tn{u) - 1)(... ®h®a® ...) (20) 
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Proof Recall the commutation relation for the Yangian: 

(1 - -^)T,{u)Uv) = T^{v)T,{u){l 



u — V u — V 

Then putting u — v = \'\ one obtains: 

A2Ti{u - h)T2(M - 2h) =T2{u- 2h)Ti(M - h)A2 

which can be rewritten as 

A2e-^^^Ti{u)e-^^^T2{u) = e'^^^T^iu - \\)e-^^^T^{u + h)A2. 

One has also to note that 

A2(Ti(n) + T2(n)) = (Ti(n) + T2{u))A2 

because 

A2{T^{u) + T2{u)) = {l-P){Ti{u)+T2{u)) = {Ti{u) + T2{u))-P{T,{u)+T2{u)) 
= (Ti(n) + T2(n)) - (Ta^ + T,iu))P = (Ti(n) + T2(n))(l - P). 

Hence we obtain 

A2((e-h^«TiH - l)(e-h^«T2H - 1)) 
= A2ie-^^-T,iu)e-^^-T2iu) - e-h^«(Ti(M) + T2(m)) + 1) 
= ie-^^-T2iu - h)Ti(M) - e-h'9„(j^^(^) + T2iu)) + 1)^2 
= (by definition) B{u,h,du)A2 (21) 

Let us denote by Aj^j+i = Id — Pj,j+i - the antisymmetrizer in i-th and i + 1-th place and 
also -Bii+i the operator i? acting on the i-th. and i + 1-th tensor components. Due to the 
formula 1^ we see that 

Recalling that the operator v4j_j+i commutes with operators not acting on the i,i + 1-th 
tensor components of [([y-'^®^ we obtain 

An{e-^^^Ti{u) -I)... {e-^^-Tn{u) - 1)(... ®a®h® ...) 

= ]^AnA,^i+^{e-^^^T^{u) - 1) . . . {e-^^^Tn{u) - 1)(... ®a®h® ...) 

= iA„(e-h^«riH - 1) . . . (e-h^«r,_i(«) - 1) 

X Ai,,;+i(e~h^«Ti(M) - l){e-^^-Ti+^{u) - 1) . . . (e~'''^"T„,(M) - 1)(... ® a ® 6 ® ...) 

= iA„(e-h^«riH - 1) . . . (e-h^«r,_i(iz) - l)fi,,+i . . . (e-h^«T„(w) - 1) 
X (g) a (g) 6(g) ...) 

= -iA„(e-h^"Ti(n) - 1) . . . (e-h^«T,_iH - l)i?,,+i . . . (e-^^^T^w) - 1) 
X (g) 6 (g) a (g) ...) 

= -A^{e-^^-Ti{u) - l){e-^^-T2iu) - 1) . . . ie~^^-TM - ® & ® a ® •••)(22) 



9 



The lemma is proved □ 

The proposition follows immediately □ 

We propose here another technical lemma which was used in the first version of the 
proof and can be useful in similar problems. 

Lemma 3 // AnX = YAn then A^XAn = nlAnX, where X,Y e A(g) Mat^'^^, A is 
an arbitrary associative algebra, Mat is the algebra of matrices and An G Mat®"' is the 
antisymmetrizer (C")*^" (C")*^"". 

Proof Recall that A^^ = A n'lT'^-i SO AfiX An — Y A^A^ — Y Anfil — nlAnX. n 

Theorem 2 Given a solution S{u) with values in [C"')®"-^V (where V is a representation 
ofU{Q[n)®^) of the equation 

An{L^{u) - 9„)(L2(n) - 9„) . . . {LM - dn)S{u) = (23) 

we take with values in V defined by 

AnS{u) = ^!{u)An{ei® ...(^en). (24) 

\1/(m) provides a solution of the equation 

QDn{u)^<{u) = 

Proof Using corollary^ we see that 

= n\An{U{u)-du){L2{u)-du)...{LM-dn)S{u) 

= An{L^{u) - d^){L,{u) -du)... {Ln{u) - d^)^AnS{u) 

= QL>„(M)^(M)A„(ei® ...®e„) 
where we have used that y4„(n!)~^ is a projector. Hence QD„(m)\1/(m) = 00 

3.2 Monodromy properties 

Let us firstly recall the fact that quantum Gaudin hamiltonians can be restricted to the 
invariant part V = (Va^ (8> ... Va^)^'^'^*'"-' invariant) tensor product of representation 

with respect to the diagonal action and this space is usually considered as a representation 
for the quantum Gaudin model. 

Hypothesis 1 The equation 

{L{u) - du)S{u) = 

for the C" V -valued function S{u) has a fundamental solution without monodromy. 
Modulo this hypothesys one obtain the following 

Theorem 3 In the case of the representation space V = {Vx^ ® ... ® Vx„)^^'^^"^ invariant) 
all the solutions of QDn{u)'^{u) = have trivial monodromy at infinity and at the points 

Z = Zi. 
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Remark 4 If we do not take the invariant part ofV\-^'^...®V\^ the monodromy at infinity 

is not zero, hut is easily controllable. 

Proof One can construct a solution to the equation 

- du){L2{u) -du)... {Ln{u) - du)S{u) = 

just taking the fundamental solution of the equation 

- du)S{u) = 

and considering S{u) = Vi® ■■■Vn-i ® S{u) where Vi are arbitrary vectors from C". Then 
one has to apply Theorem |21 to obtain solutions □ 

Remark 5 In the case of other representations like Verma modules the monodromy is 
expected to be easily controllable. In the case of sin the monodromy becomes trivial after 
the mapping GLn PGLn (recall that SL^ and PGLn are Langlands dual groups, so it 
is not accidental). 

3.3 The main correspondence 

Modulo the hypothesys [T] one obtains: 

Theorem 4 The common eigen-vectors of the Gaudin hamiltonians correspond to the 
scalar differential operators ("G-opers") hi{u)dl^ which do not have monodromy, where 
hiiu) are values of the quantum hamiltonians {—iyGnQIi{u) on the joint eigen-vector. 
Proof Following theorem IHl one just need to restrict the differential operator 

QDM = Y,{-irClQk{u)dl 

i 

to a common eigen-vector of the quantum Gaudin hamiltonians QIi{u) □ 

4 Relations 

In this section we recall several relations of our principal subject with the geometric 
Langlands correspondence of Beilinson and Drinfeld, with the separation of variables 
technique, and annonce some generalizations of the obtained results. 

4.1 Langlands correspondence 

The geometric Langlands correspondence predicts a bijection between Hitchin D-modules 
on the moduli space of vector bundles on the one side and local systems on the curve 
on the other side. Some of local systems can be realized as connections having a flag 
of invariant subbundles - such connections can be defined as differential operators called 
G-opers. On the automorphic (Hitchin) side the Hitchin D-module is just the D- module 
generated by the differential operators Hi—Ci where Hi are quantum Hitchin hamiltonians 
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and Cj are arbitrary constants parameterizing Hitchin D-modules. So there should appear 
a correspondence between some constants Ci and some differential operators (G-opers). 

The Langlands correspondence can be generalized to the case of vector bundles with 
parabolic structures, in this case instead of the Hitchin model appears the Gaudin model. 

Conjecture 2 Our quantum commuting Gaudin hamiltonians generate the same D- 
module as the hamiltonians obtained from the center of gl^ on the critical level (these 
objects are usually considered in the geometric Langlands correspondence theory). 

Let us remark that the advantage of our approach is the explicitness of proposed 
formulas in contrast to the approach based on the center of gl„ on the critical level, which 
still cannot be described exphcitly except of generators of degree 2. 

Conjecture 3 Our scalar G-oper obtained by substituting the coefficients of universal 
G-oper by some constants Gi is the same as the Langlands correspondent G-oper of the 
Hitchin D-module generated by Hi — Ci. 

4.2 Separated variables (Baxter equation) 

The relation of the G-oper to the separation of variables is as follows: in separated 
variables the wave functions of the quantum system is factorized to the product of one 
particle functions \1/ = It appears that the differential (or difference) 

equation for the wave functions 'i/'^^ (called Baxter equation) is given by the scalar G-oper, 
this is proved in Frenkel's paper for 5I2 Gaudin case and seems to be true in general. We 
hope to show it in our next publications. 

4.3 Generalizations 

The results presented here can be generalized in the following directions: firstly one can 
consider the Gaudin model for other semisimple Lie algebras g - to do this one need to 
generalize the notion of the matrix $. This can be done by the consideration the universal 
r-matrix corresponding to these Lie algebras and corresponding analogues of the Bethe 
subalgebra. 

Another generahzation is the case of XXZ and XYZ models. In these cases we expect 
the following change in the formula for the universal G-oper 

jjQXXz-naive^^-^ = TrA^{Ln{exp{v)) - exp{d^))...{Li{exp{u)) - exp{d^)). 

This formula probably needs further corrections. It is motivated by two examples - 
the Gaudin and XXX model, where we have seen the classical characteristic polynomial 
det{L{u) — A) to be "quantized" to the Universal G-oper by the rule that A is changed to 
some operator A such that the commutation relation between u and A is the quantization 
of the Poisson bracket between u and A in separated variables. Besides it is known for 
XXZ model that in separated variables the Poisson bracket {u. A} = u\. In the same vein 
one can try to generalize the formula for the universal G-oper by for XYZ model. We 
believe that our propositions about the monodromy can be generalized to these examples 
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also and hence they eventually will give a new approach to the Bethe ansatz equations in 
these models. 

Another kind of generalization concerns the quantization problem for the inte- 
grable systems with higher poles Lax operator, the so-called polynomial matrix mod- 
els |Beauville90t lAHHPj and related Hitchin models on singular curves |CT03-2j or for 
integrable systems obtained in the framework of |Chernyakov03| . The formula for the 
universal G-oper is expected to be absolutely the same. But the results about the mon- 
odromy are more problematic, possibly all the information will be encoded in the Stokes 
matrices, instead of monodromy ones, this deserves further investigations. 
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